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Abstract 

We have searched for decays that violate both lepton and baryon number using data collected 
by the Belle detector at the KEKB asymmetric e'^e" collider. No signals are found in r~ p7, 
r" PT^^, T~ Avr", and r~ — > Avr" and we set upper limits on the branching fractions 
of B{t- p-f) < 3.0 X 10-^ Bir- pir^) < 6.5 x IQ-^ B{t- Avr") < 1.3 x 10"^ and 
B{t~ Avr^) < 0.70 x 10~^ at the 90% confidence level. The former two results improve the 
previous limits by a factor of 12 and 23, respectively, while the latter two are the first searches ever 
performed. 

PACS numbers: 11.30.-j, 12.60.-1, 13.35.Dx, 14.60.Fg 
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INTRODUCTION 



While the Standard Model (SM) assumes both the baryon number [B) and lepton number 
(L) conservation, in some extensions beyond the SM such as Grand Unified Theories(GUTs), 
B and L violation is expected while their difference B — L is conserved [1] . In Ref . [2] , right- 
handed four-fermion couplings that conserve B — L were used to consider B and L violating 
r lepton, D and B meson decays. High luminosity B-factories provide an opportunity to 
look for such decays with unprecedented sensitivity. 

We report here our searches for r~ ^ p7, p7r°, Avr", and Att" decays with data samples 
of 86.7 fb~^ for jry and 153.8 fb~^ for all other modes, collected with the Belle detector at the 
KEKB asymmetric e~^e~ collider [3] . Previously, the best upper limits for the corresponding 
branching fractions were obtained by CLEO [4] based on a data sample of 4.7 fb~^: i3(r^ — >• 
p7) < 3.5 X 10"^ and B{t~ p7r°) < 1.5 x 10~^ at 90% C.L. The decay modes r~ Att" 
and r~ Att" have never been studied before. Unless otherwise stated, charge conjugate 
decays are implied throughout this paper. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cerenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter comprised of CsI(Tl) crystals (ECL) located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron fiux- 
return located outside of the coil is instrumented to detect mesons and to identify muons 
(KLM). The detector is described in detail elsewhere [5]. 

We search for r^r~ events in which one r decays into the mode studied (signal side) and 
the other r (tag side) decays into one charged particle, photons and neutrino(s). The selec- 
tion criteria are determined by examining Monte Carlo (MC) simulations for signal r-pair 
decays and background (BG) events coming from generic r-pair decays (r'^r"), qq contin- 
uum, BB, Bhabha and fifi as well as two-photon processes. The KORALB/TAUOLA [6] 
and QQ [7] generators are used for event generation, and GEANT3 [8] is used to simulate the 
Belle detector response. The two-body decays of the signal r are assumed to have a uniform 
angular distribution in the r's rest frame. All the kinematical variables are calculated in the 
laboratory frame, while those in the e+e~ center-of-mass (CM) frame are indicated by the 
superscript "CM". 
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ANALYSIS OF r" py AND pvr" 
Event selection 

The experimental signature of these events has one r lepton decaying into a proton and 
photons (signal side) and the other decaying via the 1-prong mode (tag side): 

+ nf^} + {r+ (a track)+ + nj^^^ + X(missing)} , (l) 

where a track should have transverse momentum pt > 0.1 GeV/c and polar angle —0.819 < 
cos^ < 0.906. The two highlighted tracks should have zero net charge. A proton track 
forming the signal jry/piT^ is required to have p > 1.0 GeV/c for reliable identification. A 
photon should have energy > 0.1 GeV in the same polar angle range as that for tracks. 
The number of photons on the signal side is 1 < {n^^^ + n^'^'^) < 3 for r~ p'j and 
2 < {nf'^ + n!^^^) < 4 for r~ p7r°. 

We divide the event into two hemispheres in the CM frame by the plane perpendicular 
to the thrust axis. The hemisphere containing a proton track (its identification is later 
explained), is the 'signal side', the opposite hemisphere is the 'tag side'. We require an 
invariant mass of the visible particles on the tag side of Mtag < 1-2 GeV/c^ for the py 
mode and Mtag < i^t for the pn^ one to reduce background (BG) from the e~^e~ qq 
continuum process (where q = u, d, s, c). 

For p7, we require the photon to have > 0.5 GeV. In order to reduce the number of 
fake 7's originating from a n in the qq continuum, the ratio of photon energy deposition in 
a 3 X 3 ECL square to a 5 x 5 square, E9/E25, is required to be larger than 0.93 since the 
hadronic shower of a n in the ECL is wider than a photon's electromagnetic shower. 

The 7r° from pn^ is reconstructed from a 77 pair with an invariant mass within ±5a.,^o 
((TttO = 5 — 8 MeV/c^) of the nominal 7r° mass. We impose a ir^ veto on photon(s) for Jry/pTT^ 
candidates: it should not reconstruct a 7r° meson when combined with any other photon 
whose energy exceeds 50 MeV. 

Correlations are considered among the tracks, photons and a missing particle that carries 
away undetected momentum and energy. A requirement on the total visible energy 0.5 < 
^vis I \^ < 0.92 is imposed to reject Bhabha scattering and jX^j^r production. Restrictions 
on the opening angle between p and 7, and p and 7r° reduce BG from the generic r+r~ and qq 
continuum: 0.6 < cos6'p^ < 0.96 and 0.0 < cos6'p^ < 0.95. The opening angle between the 
two tracks in the CM frame is also required to be larger than 90°. Constraints on the missing 
momentum Pmiss and polar-angle O^i^s are imposed to ensure that the missing particles are 
undetected neutrino(s) rather than photons or charged particles that escaped detection: 
Pmiss > 0.6 GeV/c for pry and Pmiss > 0.5 GeV/c for p7r°, and —0.866 < cos6'miss < 0.956 for 
both modes. A requirement on the opening angle 6'Sig'-niiss between the tag-side track and the 
missing particle helps to remove t^t~ events: cos 0^^_^:^^^ > 0.3 for Jry and cos O^^^^iss > 0-0 
for pn^. 
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FIG. 1: Pmiss vs JTT-^iss plots of T — > p^: (a) signal MC, (b) t~^t MC and (c) qq (uds) continuum 
MC. The area between the hnes is the selected region. 



An additional condition is imposed on the relation between Pmiss and the mass squared 
of a missing particle m^j^^ (see Fig.l): p^iss > -3 ml^-^^ - 1 and Pmiss > 1-2 m^j^^ - 1 for pj, 
and Pmiss > —0.8 m'^J^lss + 0-3 and Pmiss > 1-6 rri^^^^ — 1 for pn^. This cut removes 81% of the 
remaining t^t~ and 73% of continuum BG's while retaining 78% of the signal. 

Particle identification, that is very important in this measurement, is based on the re- 
sponses of subdetectors such as the ratio of the energy deposited in ECL to the momentum 
measured by CDC, the shower shape in ECL, particle range in KLM, the hit information 
from the threshold type ACC, dE/dX in CDC and time-of-fiight from TOP. We use likeli- 
hood ratios to distinguish hadron species, for instance, V{j)/tx) = Cp/{Cp + C-,,), where Ci 
is the likelihood for the detector response to the track with flavor hypothesis i. Por lepton 
identification, we use an electron probability 'P(e) and a muon probability 'P(yu) determined 
by the detector responses. 

We demand 'P(e) < 0.8 and 'P(yLt) < 0.8 for the signal-side track in order to remove 
Bhabha and fi'^fi~ processes. To identify protons, the ACC plays an important role since it 
allows us to distinguish flavors by the threshold momentum: for protons, this is p^S^ — 5 
GeV/ c in the barrel region and ~ 4 GeV/c in the forward endcap, while for pions is p^S^ > 1 
GeV/c. 

Most EG tracks are pions with a rather high momentum up to 5 GeV/c, and the mo- 
mentum of most protons ranges from 1 to 5 GeV/c. Therefore, we can identify protons and 
remove pions by requiring ACC not to fire at 1 GeV/c < p < p^S^- With these criteria, 
90% of the BG is removed, while 80% of the signal remains. 

The information from TOP and CDC is combined with that from the ACC into Cp and 
C-^ to further reduce vr BG as well as K^s in pn^. The requirement V{p/'k) > 0.8 removes 
more than 90% of the BG while retaining 87% of the signal. In the p7r° mode, background 
processes with kaons, mostly from r K*v and K* —>■ Kn^ are rejected by demanding 
V{p/K) > 0.8: 70% of the K^s are removed while 77% of the signal is retained. 
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Expected backgrounds and blind analysis 

A signal candidate is examined in the two-dimensional space of the jry / pir^ invariant 
mass, Minv, and the difference of its energy from the beam energy in the CM system, AE. 
A signal event should have Minv — i^t and AE ^ 0. The Minv and AE resolutions are eval- 
uated from the MC distributions around the peak using an asymmetric Gaussian shape to 
account for initial state radiation and ECL energy leakage for photons: a^f^^^"^ = 10.7/16.7 
MeV/c^ and cT^f Z'"" = 35.2/62.8 MeV for the p7, and = 11.3/14.9 MeV/c^ and 

^high/low _ 34 3^57 ]^ MeV for the p7r°, where the "high/low" superscript indicates the 
higher/lower side of the peak. 

To avoid any bias in extracting the result, we blind the following region: iScTMjnv 
iScTAE for p7, and the 1.68 GeV/c^ < Minv < 1-85 GeV/c^ band region for p7r° (see Fig. 2). 
Since the ±5cr region for the Jm^ mode contains significantly more events than in the Jyy 
case, we blind a wider region for the former mode. 

To estimate the expected number of BG events in the blinded region, we approximate the 
data distribution by a combination of Landau and Gaussian functions with a few parameters 
in the region 1.42 GeV/c^ < Minv < 2.10 GeV/c^ for p-f, and by an asymmetric Gaussian 
function with a few parameters in the region 1.5 GeV/c^ < Minv < 2.0 GeV/c^ for Jm^ 
as shown in Fig. 3. One finds in these areas, excluding the blinded region, 49 and 195 
data events for jry and pir^, respectively, and 51.8 ±6.1 and 178.4 ± 12.8 corresponding MC 
events (properly normalized to the luminosity of the data), after application of all cuts. The 
number of BG events expected from the BG functions in the 5a regions is 9.1 ± 1.7 and 
52.2 ± 7.3 events for p7 and pn^, respectively. The independent evaluation of the number 
of BG events in the 5a regions based on the MC simulation of the BG processes gave 6.1 
± 2.0 and 43.4 ± 6.6 events for p7 and pir^, respectively, in reasonable agreement with the 
sideband based evaluation above. The data and signal MC distributions in AE vs M^nv are 
shown for both modes in Fig. 2. 

The signal detection efficiency for the ±5a box is evaluated from MC as 9.4% and 5.8% 
for p7 and p7r°, respectively. 

Blind opening and evaluation of the branching fractions 

Data distributions after opening the blinded regions are shown in Figs. 2 and 3. We find 
16 and 70 events in the 5a regions for Jry and pvr", respectively, while the number of expected 
BG events is 9.1 ± 1.7 and 52.2 ± 7.3 as evaluated above using the data sidebands. The 
differences between the number of observed data events and the BG expectations, 6.9 ± 4.3 
for pj and 17.8 ± 11.1 for pn^, are not statistically significant. 

To extract the number of signal events, we applied an unbinned extended maximum 
likelihood method (UEML), which is more sensitive than the binned one since it uses the 
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FIG. 2: Ai? - Minv plot after opening the blind areas (inside the dotted hnes) for (a) p'j and (b) 
pTT^ decay modes. The data are shown by bold dots and the signal MC by the boxes. 



complete information about events. In this method the likelihood function is defined as 

-{s+b) N 

C{s,b) = ——l[{sS. + bB,), (2) 
1=1 

where N is the number of the observed events, s and b are the free parameters corresponding 
to the number of signal and background events, respectively, and Si and Bi are the values 
of the probability density functions of the signal and BG for the z-th event. 5*^ is given by 
the signal MC distributions, and Bi is the BG function obtained above and normalized to 
unity. 

The maximum likelihood fit for the ±5cr region yields Sq = 0.16 and b^ = 15.84 for pj 
and So = 3.09 and bo = 66.91 for pir^, respectively. Following Ref. [9], the upper limit at 
the 90% confidence level (C.L.) is obtained by means of toy MC, as described below. For 
every assumed expected signal yield s, 10,000 samples are generated, for each of which the 
number of signal and BG events is determined by Poisson statistics with the mean values s 
and bo, respectively. We then assign the AE and Mjnv values to these events according to 
their distributions; a UEML fit is performed for every prepared sample to extract the signal 
yield (sq); the confidence level for an assumed s is defined as the fraction of the samples 
whose §0 exceeds s. This procedure is repeated until we find the value of sgo that gives a 
90% chance of sq being larger than s. The resulting values are sgo = 3.7 events for p7 and 
Sgo = 9.8 events for pir^. 

The upper limit on the branching fraction B at the 90% C.L. is then calculated as 
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FIG. 3: The AE distributions for (a) p7 and (b) pvr". In both (a) and (b) the upper, middle and 
lower figures correspond to the indicated Minv regions: the middle one is the blind area while the 
upper and lower correspond to the sidebands. All three regions have the same width. The closed 
circles with error bars are the data observed after opening the signal box, the histogram with error 
rectangles shows the results of MC simulation for the BG processes (generic t~^t~ and continuum 
qq). The thick solid curve shows the contribution of the BG parameterized in the sideband regions. 
The expected signal (^B = 10 x 10~^) is shown by the thin solid line. The data in the middle figures 
are blinded until the last stage of the analysis. 



where the number of produced r-pairs is A^^,- = 78.9 x 10^ and 140.0 x 10^ for and pyr", 
respectively, and the detection efficiencies are e = 9.4% and 5.8%. The resulting upper limits 
are B{t- p-f) < 2.5 x 10"^ and B{t- p7r°) < 6.1 x 10"^ at 90% C.L. 

Systematic uncertainties related to the detector sensitivity 2eNrT in the denominator of 
Eq. (3) are evaluated to be 5.9% for p-f and 8.7% for pn^. The individual contributions to the 
uncertainties for Jry (pir^) are 2.0%(2.0%) from tracking efficiency, 3.0%(6.0%) from photon 
reconstruction efficiency, 2.0%(4.0%) from selection criteria, 3.0%(3.0%) from trigger effi- 
ciency, 2.5%(3.0%) from proton identification, 0.2%(0.3%) from MC statistics, 1.4%(1.4%) 
from luminosity evaluation, and 0.03%(0.03%) from the r"^r~ cross-section. 

The systematic uncertainty in sgo is estimated by varying the parameters of the BG 
functions by ±lcr. This gives uncertainties of ±0.77 events for p7 and ±0.75 events for 
pn^. The upper limits on the branching fractions taking into account all the systematic 
uncertainties are calculated to be 



B{t~ p-i) < 3.0 X 10~^ 
B{t' ^p7r°) < 6.5 X 10-^ 



(4) 
(5) 



10 



at 90% C.L. with 86.7 fb~^ and 153.8 fb~^ of data, respectively. 

The p7r° mode should be discussed briefly since some excess of events in the signal area 
can be seen in the AE distribution of Fig. 3. From the s dependence of the likelihood 
function, a 68% confidence interval for sq is —1.6 < sq < 8.8; in other words, Sq is consistent 
with within 1 a. If we use toy MC, a probability to have sq > 3.09 when the signal yield 
5 = is 19%. We conclude that the observed excess can be due a statistical fluctuation. 

ANALYSIS OF r" Att" AND Att" 
Event Selection 

The experimental signature of these events has one r lepton (signal side) decaying to Avr, 
A to pvr and the other (tag side) decaying via a 1-prong mode: 

{r" ^ (p7r+) + TT"} + |r+ ^ (a track)+ + (n™*^ > 0) + X(missing)|. 

We consider here the B — L conserving decay modes: r~ Avr^ with A pvT"*" and 
r+ — s> Att''" with A pn^ . We also consider the B — L violating decay modes: Avr" 
with A p-K" and Att"*" with A — p7r+. The experimental signature for these modes 

is: 

{r" ^ (pvr") + TT"} + |r+ ^ (a track)+ + (n!^"^*^ > 0) + X(missing)|. 

We denote the pion from r — Avr as tti and the pion from A ^ pvr as 7r2. We can distinguish 
between the B — L conserving and violating modes by the charge of these pions: the B — L 
conserving decay modes have an opposite sign combination on the tti and charges, while 
the B — L violating modes have a same sign combination. As in the case of — > jrf and 
p7r° modes, tracks and photons should have > 0.1 GeV/c and > Q.l GeV, respectively, 
with a polar angle satisfying —0.866 < cos^ < 0.956. 

We first demand that the four tracks have a net zero charge. The magnitude of the thrust 
is required to be larger than 0.9 to suppress the qq continuum background. The event should 
have a 1-3 prong configuration relative to the plane perpendicular to the thrust axis. We 
select A candidates via the pvT"*" decay channel based on the angular difference between the 
A flight direction and the direction pointing from the interaction point to the decay vertex 
(see, Ref.[10] for more details). The proton from the A decay is identified by demanding 
V{p/7t) > 0.6. In order to avoid fake A candidates in which two tracks in the signal side are 
an e~e~^ pair from a photon conversion, the electron veto is imposed on the three tracks in 
the signal side. The reconstructed A candidate mass should be within ±5 MeV/c^ of the 
nominal A mass and > 1.75 GeV/c is required to reduce contributions from the generic 
T^T~ and qq continuum as shown in Fig. 4. 

As in the p-f and pir^ cases, the following criteria are imposed: 5.29 < E^^ < 10.5 GeV to 
reject Bhabha, two-photon and fj,~^fj.~ reactions; Pmiss > 0.4 GeV/c and —0.866 < cos 6*111188 < 
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Mass (GeV/c^) p^cM (GeV/c) 



FIG. 4: Reconstructed A candidate mass (left) and momentum (right) in the CM frame. The 
signal MC distributions are indicated by the filled histogram, all BG's including r~^r~ and qq by 
the open histogram, and closed circles are data. While the signal MC is normalized arbitrarily, the 
data and MC are normalized to the same luminosity. The selected area is indicated by the lines 
with arrows. 

0.956 to ensure that a missing particle is a neutrino(s); cos 0^^_^^^ > to include the missing 
particle in the tag side; n^'^ < 1 and n^^'^ < 2 to suppress the continuum background. 

Both the proton veto V{p/n) < 0.6 and kaon veto V{K/n) < 0.6 are applied to tti and 
the tag-side track. 

The correlation between the missing momentum Pmiss and mass-squared m^^gg is consid- 
ered to further suppress BG's from generic t^t~ and continuum BG: Pmiss > 1.5 xm^j^^ — 1.0. 

Signal resolutions and blind analysis 

The Minv and AE resolutions are evaluated from MC: o^f^j^^^"" = 4.6/4.0 MeV/c^ and 
Cae''''°™ = 22/29 MeV. In Avr" (Avr"), there is no M^^^ tail due to energy leakage from ECL 
because there are no photons in the final state of this mode. 

We blind a region over ±5aMi^^ and —0.5 < AE < 0.5 GeV. Fig. 5 shows scatter-plots 
for data and MC samples over ±15(T in the Mjnv — AE plane: the number of data and MC 
events outside the blinded region (bounded by the vertical dotted line in Fig. 5 (a) and (b)) 
are 11 and 13.2 ± 3.5 events, respectively. Good agreement is observed. The surviving BG 
events are due to generic t~^t~ decays (about 1/2) and uds continuum (about 1/2). The 
former events are dominated by the r — > 01(1260)//,- decays, in which three charged pions 
from the ai(1260) decay form a fake A candidate. The continuum BG events have one true 
A which forms a signal candidate together with another track. 
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FIG. 5: Scatter-plot of data and MC events in the Mjnv ^ Ai? plane: (a) and (b) correspond to the 
±15(7 area for the B — L conserving and violating modes, respectively, while (c) and (d) show the 
region 1.5 < M\av < 2.0 GeV/c^ and |AE'| < 0.5 GeV for the same modes. The areas of (a) and 
(b) are also shown as solid rectangles in (c) and (d). The areas inside the dotted lines in (c) and 
(d) denote the blind signal regions. The 90% elliptical region shown by a solid curve in (a) and 
(b) is used for evaluating the signal yield. In (a) and (b) the vertical dotted lines denote the blind 
regions similar to those in (c) and (d) ; the regions inside the horizontal dotted lines and outside the 
vertical dotted lines are sidebands used to estimate the expected BG in the elliptical region. Closed 
circles correspond to the data (154 fb~^), open squares are generic MC t~^t~ events (equivalent 
luminosity of 297 fb~^) and open triangles are MC uds continuum events (equivalent luminosity 
of 104 fb"'^). Filled boxes show the MC signal distribution with arbitrary normalization. 



13 



Blind opening and evaluation of the branching fraction 



Since there are fewer remaining events compared to Jry and p7r°, we apply the Frequentist 
approach using the Feldman & Cousins method [11] rather than the maximum hkehhood 
method. We take an elhptical region that contains 90% of MC signal events passing all 
cuts as a signal region, as shown in Fig. 5 (a) and (b). It results in the signal detection 
efficiency e = 11.8% for the B — L conserving and e = 11.7% for the B — L violating modes, 
respectively. 

From Fig. 5 (c) and (d) we assume the BG distribution to be flat along the Minv axis, 
and then obtain the expected BG in the ellipse as 1.7 ±0.8 events for each of the two modes, 
using sideband regions, inside the horizontal dotted lines and outside the vertical dotted 
lines, as shown in the Fig. 5 (a) and (b). We open the blinded region and flnd only one data 
event in the ellipse for the B — L conserving mode and no data events for the B — L violating 
mode (see Fig. 5 (a) and (b), respectively). The upper limits on the signal yields at 90% 
C.L. are obtained by the Feldman- Cousins method as sgo = 2.8 and sgo = 1.2, respectively. 
The upper limits on the branching fraction are then calculated as 

where N^r = 140.0 x 10^ and 13{A pn) = 0.639 [12]. The resulting values are B{t- 
Att") < 1.3 X 10~^ and B{t~ Avr") < 0.58 x 10"^. 

Among systematic uncertainties on the detection sensitivity 5*0 = 2eNT-TB{A — » pn), A 
selection and the proton identiflcation in A decay contribute 6.0% and 3.0%, respectively; 
B{A pn) has an uncertainty of 0.8% [12]; the trigger efficiency (0.5%), tracking (4.2%), 
selection criteria (4.0%), MC statistics (0.7%), luminosity (1.4%) and the r+r" cross-section 
(0.03%) are also considered. All these uncertainties are added in quadrature to 9.1% in total. 

The upper limits on the branching fractions at the 90% C.L. including systematic errors 
are then calculated by the POLE program [13]. The resulting branching fractions are 

B{t^ Att") < 1.3 X 10"^ 
B{r^ An-) < 0.70 x 10~^ 

at the 90% C.L. 



RESULTS 

We obtain the following preliminary upper limits on the branching fractions: B{t~ 
P7) < 3.0 X 10-^ B{r~ pn^) < 6.5 x 10"^ i3(r~ An~) < 1.3 x 10"^ and B{t- 
An-) < 0.70 X lO"'^ at the 90% confidence level. 

For the latter two modes this is the first search ever performed. 
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The resulting upper limits on the branching fraction for r~ jry and r~ pir^ improve 
upon the previous measurements by a factor of 12 and 23, respectively. This large improve- 
ment is mostly due to the powerful proton identification ability of the Belle detector that 
removes spurious combinatorial BG's as well as higher statistics compared to the previous 
experiment. 
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